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ABSTRACT 

We report the results of spectrophotometric observations of the massive star MN18 
revealed via discovery of a bipolar nebula around it with the Spitzer Space Telescope. 
Using the optical spectrum obtained with the Southern African Large Telescope, we 
classify this star as B1 la. The evolved status of MN18 is supported by the detection of 
nitrogen overabundance in the nebula, which implies that it is composed of processed 
material ejected by the star. We analysed the spectrum of MN18 by using the code 
CMFGEN, obtaining a stellar effective temperature of Ri21 kK. The star is highly red¬ 
dened, E{B — U)«2 mag. Adopting an absolute visual magnitude of My=—6.8±0.5 
(typical of B1 supergiants), MN18 has a luminosity of log L/ Lq«5.42±0.30, a mass- 
loss rate of Ri(2.8 — 4.5) x 10“^ M 0 yr“^, and resides at a distance of kpc. 

We discuss the origin of the nebula around MN18 and compare it with similar nebulae 
produced by other blue supergiants in the Galaxy (Sher25, HD 168625, [SBW2007] 1) 
and the Large Magellanic Cloud (Sk—69°202). The nitrogen abundances in these neb¬ 
ulae imply that blue supergiants can produce them from the main sequence stage up 
to the pre-supernova stage. We also present a AT-band spectrum of the candidate lu¬ 
minous blue variable MN56 (encircled by a ring-like nebula) and report the discovery 
of an OB star at ~17 arcsec from MN18. The possible membership of MN18 and the 
OB star of the star cluster LyngaS is discussed. 

Key words: circumstellar matter - stars: emission-line. Be - stars: individual: 
[GKF2010] MN18, [GKF2010] MN56, [GKF2010] MN109 ~ supergiants - stars: winds, 
outflows open clusters and associations: individual: LyngaS. 


* Based on observations made with the Southern African 
Large Telescope, programs 2011-3-RSA_OTH-002 and 
2013-1-RSA_OTH-014, and the 2.2m MPG telescope at La 
Silla (ESO), program ID P092.A-9020. 
t E-mail: vgvaram@mx.iki.rssi.ru 


1 INTRODUCTION 

Massive stars change their mass-loss rate, and the velocity, 
symmetry and composition of their copious winds in the 
course of their evolution. The mass loss along with stellar 
rotation and magnetic fields, and effects of binarity (various 
modes of mass transfer between binary components, merg¬ 
ers of the components, common-envelope ejection) result in 
the formation of different types of circumstellar nebulae, of 
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which the hourglass-like type is most spectacular. The first 
and the best known example of such a nebula associated 
with a massive star is that created by the Large Magellanic 
Cloud B-type supergiant Sk—69°202 - the progenitor of the 
SN 1987A (Plait et al. 1995; Burrows et al. 1995). Later, a 
number of Galactic analogs of this nebula were discovered 
around the already known blue supergiants - the candidate 
luminous blue variables (cLBVs) Sher25 (Brandner et al. 
1997a), HD 168625 (Smith 2007) and MWC349A (Gvara¬ 
madze & Menten 2012), and the newly identified blue su¬ 
pergiant [SBW2007] 1 (Smith, Bally & Walawender 2007). 

There is still no complete understanding on how and in 
which evolutionary phases massive stars form their circum- 
stellar nebulae. The nebulae can definitely originate during 
the late and final (pre-supernova) phases, as evidenced by 
the detection of numerous compact shells of processed ma¬ 
terial around Wolf-Rayet stars (e.g. Chu, Treffers & Kwit- 
ter 1983; Esteban et al. 1992; Dopita et al. 1994; Marston 
1995; Stock, Barlow & Wesson 2011) and the three-ring neb¬ 
ula associated with the SN 1987A. Under proper conditions, 
compact shells could form around red supergiants (Mackey 
et al. 2014), which unlike the shells around Wolf-Rayet stars 
could exist until their central stars explode as supernovae. 
The majority of candidate and bona fide LBVs are also as¬ 
sociated with compact circular or bipolar nebulae (Nota et 
al. 1995; Glark, Larionov & Arkharov 2005; Kniazev, Gvara¬ 
madze & Berdnikov 2015) and there is evidence that some of 
these stars might be the immediate precursors of supernovae 
(e.g. Kotak & Vink 2006; Groh, Meynet & Ekstrom 2013; 
Justham, Podsiadlowski & Vink 2014). On the other hand, 
the analysis of the heavy element abundances in circumstel- 
lar nebulae around several LBVs by Lamers et al. (2001) led 
them to suggest that these nebulae were created at the end 
of the main sequence. Similarly, it was argued that enhanced 
mass loss at the beginning of the blue supergiant phase was 
responsible for the origin of the bipolar nebula around the 
cLBV Sher25 (Smartt et al. 2002; Hendry et al. 2008) and 
a multi-shell nebula associated with the late O-type super¬ 
giant TYG 3159-6-1 (Gvaramadze et al. 2014a). Which of 
the evolutionary phases are more prone to produce circum- 
stellar nebulae in general and the bipolar ones in particular 
still needs to be determined. Glarifying this point is impor¬ 
tant for estimating the production rate of different types of 
nebulae and for the comprehension of the mechanisms re¬ 
sponsible for their origin. 

A large number of models have been proposed to ex¬ 
plain the formation of bipolar nebulae. The vast majority 
of them assume that the spherically symmetric stellar wind 
is collimated by a disk-like (equatorial) circumstellar envi¬ 
ronment, whose origin and geometry could be attributed to 
either an asymmetric mass loss during the preceding stellar 
evolutionary phase (e.g. Blondin & Lundqvist 1993; Mar¬ 
tin & Arnett 1995; Frank, Balick & Davidson 1995; Langer, 
Garfa-Segura & Mac Low 1999) or to various binary inter¬ 
action processes (e.g. Morris 1981; Fabian & Hansen 1979; 
Mastrodemos & Morris 1999; Morris & Podsiadlowski 2009). 
In some models, the stellar wind is considered to be intrin¬ 
sically bipolar, e.g. because of nearly critical stellar rotation 
(Owocki & Gayley 1997; Maeder & Desjacques 2001; Ghita 
et al. 2008). It is quite likely that several mechanisms for 
shaping bipolar nebulae could be at work. The detection of 
new examples of bipolar outflows from massive stars might 


be essential for the identification of relevant mechanisms for 
their origin as well as for a better understanding at which 
stellar evolutionary phase(s) these mechanisms start to op¬ 
erate. 

In our search for compact mid-infrared nebulae in the 
archival data of the Spitzer Space Telescope (for the moti¬ 
vation and some results of this search, see Gvaramadze et 
al. 2009, 2010a,b, 2014b; Gvaramadze, Kniazev & Fabrika 
2010c, hereafter Paper I), we detected two striking examples 
of hourglass-like nebulae, which were presented for the first 
time and named MN13 and MNl^ in Paper I (see tablet 
and figs3 and 10 there). MN18 and its central star (in the 
following, we will use the name MN18 only for the star) 
are the subject of this paper and presented in Section^ 
Section[3] describes our spectroscopic and photometric ob¬ 
servations of MN18. The spectral classification of MN18 is 
given in Section|4l where we also determine the fundamen¬ 
tal parameters of this star using the radiative transfer code 
CMFGEN and estimate its luminosity. Spectroscopy of the 
nebula is discussed in Section[5l In Section[6l we discuss the 
origin of bipolar nebulae around evolved massive stars and 
compare the nebula around MN18 with other similar bipolar 
nebulae. The possible parent cluster of MN18 is discussed 
in Section|5]as well. 


2 MN18 AND ITS BIPOLAR NEBULA 

The blue supergiant MN18 is one of many dozens of evolved 
massive stars revealed through the detection of their circum¬ 
stellar nebulae with the Spitzer Space Telescope (see Table 1 
in Paper I). The bipolar nebula around MN 18 was discovered 
in archival 24 /rm data obtained with the Multiband Imag¬ 
ing Photometer for Spitzer (MIPS; Rieke et al. 2004) within 
the framework of the 24 and 70 Micron Survey of the Inner 
Galactic Disk with MIPS (Garey et al. 2009). This nebula 
appears as two lobes protruding for about 70 arcsec in the 
northwest and southeast directions from a bright waist cen¬ 
tred on MN18. At a distance of d=5.6 kpc (see Section|4]), 
the linear extent of the lobes is Ril.9 pc. 

Fig.[T]shows the nebula in three different intensity scales 
to highlight some details of its structure. The low intensity 
image of the waist of the nebula revealed two bright spots on 
each side of MN18 (the star itself is not visible at 24 pm). We 
interpret these spots as cross-sections of a toroidal equatorial 
ring around the star visible nearly edge-on. The northwest 
lobe appears to be opened along the polar axis, while in the 
southeast lobe one can discern two incomplete ellipse-like fil¬ 
aments oriented plane parallel to the equatorial ring. These 
filaments are reminiscent of a polar ring associated with the 
cLBV HD 168625 (see fig. Id in Smith 2007 and fig. 2n in 
Paper I), but both are located on the same side of the torus. 
Assuming rotational symmetry for the nebula, the aspect ra¬ 
tios of the filaments suggest that the nebula is inclined with 
respect to the plane of sky by «30°. Gorrespondingly, the 
inclination angle of the stellar rotational axis to our line of 
sight, i, is «60°. Fig.[T] also shows a bright feature attached 


^ In the Set of Identifications, Measurements and Bibliography 
for Astronomical Data (SIMBAD) data base these shells are 
named [GKF2010] MN13 and [GKF2010] MN18. 
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Figure 1. Spitzer MIPS 24 fira image of the bipolar nebula around MN18 at three intensity scales, highlighting various details of its 
structure. The coordinates are in units of RA (J2000) and Dec. (J2000) on the horizontal and vertical scales, respectively. At an adopted 
distance of 5.6 kpc, 1 arcmin corresponds to Ril.6 pc. 
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Figure 2. From left to right, and from top to bottom: Spitzer MIPS 24 fim, IRAC 8, 5.8, 4.5 and 3.6^m, and SHS Ha+[N ll] images of 
the equatorial ring around MN18. MN18 and several stars around it were proposed by Lynga (1964a) to be members of the same star 
cluster (see Section r6.4l for more details). 


to the southwestern edge of the southeast lobe. A possible 
origin of this feature is discussed in Section[6]4] 

MN18 and its equatorial ring are clearly visible in all 
(3.6, 4.5, 5.8 and 8.0 /rm) images (see Fig.[2| obtained with 
the Spitzer Infrared Array Camera (IRAC; Fazio et al. 2004) 
within the Galactic Legacy Infrared Mid-Plane Survey Ex¬ 
traordinaire (CLIMPSE; Benjamin et al. 2003). The north¬ 
east half of the ring is «20 per cent brighter than the south¬ 
west one. The semimajor and semiminor axes of the nebula 
are «11 and 4.5 arcsec, respectively. For d=5.6 kpc, the lin¬ 
ear radius of the waist is «0.29 pc. The 8 pm image also 
shows a “horn” at the northeast edge of the equatorial ring, 
which might correspond to the beginning of the northwest 
lobe. Similar horns could also be seen in the IRAC images 


of the equatorial rings of the cLBVs Sher25 and MN56 (see, 
respectively, fig. 21 and fig. 4 in Paper I). 

The IRAC images of the waist of the bipolar neb¬ 
ula were presented for the first time by Kwok et al. 
(2008), who classified it as a planetary nebula (PN) and 
named PNG 321.0—00.Kwok et al. (2008) also noted that 
PNG 321.0—00.7 has an optical counterpart in the Super- 
GOSMOS H-alpha Survey (SHS; Parker et al. 2005), but 


^ Anderson et al. (2012) used mid- and far-infrared colours of 
PNG 321.0—00.7 to suggest that this object is likely a PN, while 
Cohen et al. (2011) confused it with another PN candidate 
from Kwok et al. (2008), PNG 306.4-1-00.2, and proposed that 
PNG321.0—00.7 might be a reflection nebula (see their table!). 
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Table 1. Details of MN18. 


Spectral type 

B1 la 

RA(J2000) 

15‘'16"'41!01 

Dec.(J2000) 

-58°22'26'.'0 

1 

321?0239 

b 

-0?6992 

B (mag) 

15.34±0.01 

V (mag) 

13.40±0.01 

Ic (mag) 

10.91±0.01 

J (mag) 

9.22±0.02 

H (mag) 

8.72±0.02 

Ks (mag) 

8.38±0.02 

[3.6] (mag) 

8.02±0.04 

[4.5] (mag) 

7.99±0.06 

[5.8] (mag) 

7.88±0.05 

[8.0] (mag) 

8.01±0.14 


did not show its image. Fig. 10 in Paper I presented the 
SHS Ha+[Nii] image of the waist for the first time (see 
also Fig.[2|. 

In Paper I, we provisionally classified MN18 as an early 
B supergiant, using a spectrum obtained with the 1.9-m 
telescope of the South African Astronomical Observatory 
(SAAO) (see Section lSTTI) . Our subsequent spectroscopy with 
the Southern African Large Telescope (SALT) allowed us to 
refine this classification to Blla (see Section l4Tl) . 

The details of MN18 are summarized in Table[T] The 
BVIc magnitudes are the mean values derived from our 
hve CCD measurements listed in Table[3] The coordinates 
and the JHKa photometry are taken from the Two-Micron 
All Sky Survey (2MASS) All-Sky Catalog of Point Sources 
(Cutri et al. 2003). The IRAC photometry is from the 
GLIMPSE Source Catalog (I -I- II + 3D) (Spitzer Science 
Center 2009). 

MN18 is one of the 25 stars with P < 15 mag, which 
were suggested by Lynga (1964a) to be members of the 
open star cluster Ly 3, called Lynga 3 in the SIMBAE0 and 
WEBDA|3 (Mermilliod 1995) data bases. MN18 has the se¬ 
quence number 22 in table V of Lynga (1964a), which gives 
photographic B and V magnitudes of 15.6 and 13.3, respec¬ 
tively, for this star. The relationship between MN18 and 
Lynga 3 is further discussed in Section[6]4] 


3 MN18: OBSERVATIONS 

3.1 Spectroscopic observations 

The Erst spectrum of MN18 was obtained with the SAAO 
1.9-m telescope on 2009 May 27. The observations were per¬ 
formed with the Cassegrain spectrograph using a slit of 3 
arcmin x 1.8 arcsec and grating with 300 lines mm~^. The 
slit was oriented at a position angle (PA) of PA=80°, mea¬ 
sured from north to east. This spectral setup covered a wave¬ 
length range of «4200-8100 A with a reciprocal dispersion 
of «2.3 Apixel”^ and a spectral resolution full width at half¬ 
maximum (FWHM) of «7A. Three exposures of 900 s were 
taken with a seeing of r; 2 arcsec. 

® http://simbad.u-strasbg.fr/simbad/ 
http://webda.physics.muni.cz/ 


The obtained spectrum (see fig. 9 of Paper I) revealed 
numerous H and He l lines in absorption and prominent neb¬ 
ular emissions - the Ha line flanked by the [N ll] AA6548, 
6584 lines. Unfortunately, the stellar and the nebular spec¬ 
tra were too weak to allow us to obtain a reliable spectral 
classification of MN18 and to derive any useful information 
on the nebula. 

To classify MN18 and to study the nebula, we obtained 
a new spectrum (we call it low-resolution hereafter) with the 
SALT (Buckley, Swart & Meiring 2006; O’Donoghue et al. 
2006) on 2012 March 2, using the Robert Stobie Spectro¬ 
graph (RSS; Burgh et al. 2003; Kobulnicky et al. 2003) in 
the long-slit mode with a slit of 8 arcminxl.5 arcsec. The 
slit was oriented along the major axis of the optical neb¬ 
ula, i.e. at PA=72°. The PG900 grating was used to cover 
the spectral range of 4200—7300 A with a final reciprocal 
dispersion of 0.97Apixel“^ (FWHMr;5.50 A). A 300 s spec¬ 
trum was taken with 2x2 binning. The RSS uses a mosaic 
of three 2048x4096 CCDs and the hnal spatial scale for this 
observation was Cf.'25pix“^. The seeing during the observa¬ 
tion was «1.3 arcsec. A Xe lamp arc spectrum was taken 
immediately after the science frame. A spectrophotometric 
standard star was observed during twilight time for relative 
flux calibration. 

To study the velocity distribution along the optical neb¬ 
ula, two additional higher resolution long-slit spectra (high- 
resolution spectra hereafter) were obtained with the RSS us¬ 
ing the grating PG2300 on 2013 April 26 and 27. A slit with 
a width of 1.25 arcsec was used in both observations with 
the same PA as earlier. A spectral range of 6035—6900 A was 
covered with a final reciprocal dispersion of 0.27Apixel“^ 
and a FWHM spectral resolution of 1.33 A. One 900 s spec¬ 
trum was taken with a 2x4 binning each night with a spa¬ 
tial scale of 0/51 pix~^ along the slit. The seeing during the 
observations was 1.5 arcsec. A Ne lamp arc spectrum and 
a set of Quartz Tungsten Halogen flats were taken imme¬ 
diately after the science frames, and a spectrophotometric 
standard star was observed during twilight time for relative 
flux calibration. 

The primary reduction of the SALT data was carried 
out with the SALT science pipeline (Crawford et al. 2010). 
After that, the bias and gain corrected and mosaiced long- 
slit data were reduced in the way described in Kniazev et al. 
(2008). Absolute flux calibration is not feasible with SALT 
because the unfilled entrance pupil of the telescope moves 
during the observation. However, a relative flux correction 
to recover the spectral shape was done using the observed 
spectrophotometric standard. 

For spectral modelling and to derive the rotational ve¬ 
locity, we obtained high-resolution spectra of MN18 with the 
2.2-m Max Planck Gesellschaft telescope at La Silla (Eu¬ 
ropean Southern Observatory; ESO). The data were taken 
on 2014 March 8 and 9 with the Fibre-fed Extended Range 
Optical Spectrograph (FEROS; Kaufer et al. 1999) - a high- 
resolution Echelle spectrograph with an aperture of 2 arcsec, 
which provides a resolution of 0.038 A pixeP^. A 2048x4096 
GCD was used to record 39 orders covering the spectral 
range of 3500-9500 A. Nine 1800 s exposures were obtained, 
of which three were taken during the first night. The seeing 
during the observation was below 1.5 arcsec in both nights. 
Wavelength calibration frames were obtained with ThArNe 
lamps. A spectrophotometric standard star was observed di- 
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Table 2. Journal of the observations. 


Spectrograph 

Date 

Exposure 

(min) 

Spectral scale 
(A pixel” ^) 

Spatial scale 
(arcsec pixel”^) 

Slit 

(arcsec) 

Spectral range 

(A) 

Cassegrain (SAAO 1.9-m) 

2009 May 27“ 

3x15 

2.3 

1.4 

1.8x180 

4200-8100 

RSS (SALT) 

2012 March 2 

5 

0.97 

0.25 

1.5x480 

4200-7300 

RSS (SALT) 

2013 April 26 

15 

0.27 

0.51 

1.25x480 

6035-6900 

RSS (SALT) 

2013 April 27 

15 

0.27 

0.51 

1.25x480 

3500-9500 

FEROS (MPG 2.2-m) 

2014 March 8 

3x30 

0.038 

- 

- 

3815-5440 

FEROS (MPG 2.2-m) 

2014 March 9 

6x30 

0.038 

- 

- 

3815-5440 


Note. “Presented for the first time in Paper I. 


rectly before the target observation in order to perform a 
flux calibration. 

The reduction of the FEROS data, including flatfield 
correction, bias and background subtraction, order tracing, 
wavelength and flux calibration, was performed manually 
with the FEROS Data Reduction Software based on ESO- 
MIDAS. 

The log of our spectroscopic observations is listed in 

TableEl 

3.2 Photometry 

To detect possible photometric variability of MN18, we de¬ 
termined its B,V and Ic magnitudes on CCD frames ob¬ 
tained with the 76-cm telescope of the SAAO during our 
five observing runs in 2009-2014. We used an SBIG ST- 
lOXME CCD camera equipped with BVB filters of the 
Kron-Cousins system (see e.g. Berdnikov et al. 2012). We 
also calibrated the 1.9-m telescope spectrum and synthe¬ 
sized its V magnitude in the way described in Kniazev et al. 
(2005). Additionally, the V magnitude of MN18 was mea¬ 
sured on V band images taken before each SALT spectro¬ 
scopic observation. For this we used secondary photometric 
standards calibrated with the 76-cm telescope data. Using 
the same data, we also re-calibrated the photographic B and 
V magnitudes of MN18 from Lynga (1964a). The resulting 
photometry is presented in Tabled 

From Table[3]it follows that MN18 did not experience 
major changes in its brightness: during the last five years 
its B, V and B magnitudes changed by <0.1 mag, while its 
B — V and B — B colours remained almost constant with 
mean values of 1.94±0.01 and 4.43T0.01 mag, respectively. 


4 MN18: SPECTRAL ANALYSIS AND 

STELLAR PARAMETERS 

4.1 Classification of MN18 and changes in the Ha 
line 

The normalized low- and high-resolution SALT spectra of 
MN18 are presented in Fig.[3l with the principal lines and 
most prominent diffuse interstellar bands (DIBs) indicated. 
All wavelengths are given in air. Equivalent widths (EWs), 
FWHMs and heliocentric radial velocities (RVs) of the main 
lines in the low-resolution spectrum (measured applying the 
MIDAS programs; see Kniazev et al. 2004 for details) are 
summarized in Table|4l 

Fig.[3]shows that the spectrum of MN18 is dominated by 


Table 3. Photometry of MN 18. 


Date 

B 

V 

B 

1963 April 4-17(“) 

15.46±0.25 

13.28±0.21 

- 

2009 May 27('’) 

- 

13.35±0.03 

- 

2010 January 22^*^^ 

15.36±0.01 

13.42±0.01 

10.94±0.01 

2012 March 2(‘*) 

- 

13.33±0.04 

- 

2012 May 6 (“) 

15.30±0.01 

13.36±0.01 

10.87±0.01 

2013 January 

15.34±0.01 

13.41±0.01 

10.90±0.01 

2013 April 26('^) 

- 

13.44±0.04 

- 

2014 January 18(“) 

15.33±0.02 

13.43±0.01 

10.92±0.01 

2014 April 8 ^) 

15.38±0.02 

13.41±0.02 

10.94±0.02 


Note: “Lynga (1964a); ^1.9-m telescope; “76-cm telescope; 
‘^SALT. 


Table 4. EWs, FWHMs and RVs of main lines in the low- 
resolution spectrum of MN18. 


Ao(A) Ion 

EW(A) 

(A) 

FWHM(A) 

(A) 

RV 

(kms”^) 

4340 H 7 

1.61±0.04 

8.68±0.41 

-55±9 

4388 He i 

0.26±0.02 

8.46±0.67 

- 22±8 

4471 He I 

0.94±0.04 

7.03±0.36 

-51±8 

4481 Mg II 

0.18±0.02 

5.33±0.90 

2±8 

4553 Si III 

0.56±0.02 

6.34±0.23 

-63±8 

4713 He I 

0.18±0.01 

4.85±0.21 

- 11±8 

4861 H/3 

1.61±0.05 

9.19±0.31 

-48±8 

4922 He I 

0.59±0.01 

6.40±0.17 

-37±7 

5045 N II 

0.31±0.01 

8.85±0.39 

-44±7 

5667 N II 

0.53±0.02 

8.13±0.18 

- 68±6 

5680 N II 

0.82±0.02 

8.13±0.18 

-104±6 

5740 Si III 

0.67±0.01 

7.30±0.17 

-34±6 

5876 He i 

0.90±0.02 

6.91±0.17 

-78±6 

6563 Ha(“) 

0 . 20 ± 0.01 

6.70±0.21 

-221±5 

6678 He i 

1.15±0.02 

6.51±0.15 

-69±6 

7065 He i 

1.13±0.04 

9.41±0.41 

-90±5 


Note: “The RV of this line is noticeably affected by the P Cygni 
absorption. 

absorption lines of H and Hel. In the high-resolution spec¬ 
trum the Ha line is in emission and shows a P Cygni profile. 
Other prominent absorption lines are those of N ll AA5001-5, 
5045, 5667, 5680, 6380, 6482, On AA4350, 4415-17, 4640-50, 
6721, Mg I A4471, Sim AA4553-68-75, 4650 5740, and Cii 
AA6578, 6583. No Hen lines are present in the spectrum, 
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Figure 3. Normalized low-resolution spectrum of MN18 obtained with the SALT on 2012 March 2 with the principal lines and most 
prominent DIBs indicated. The lower panel also shows the high-resolution spectrum obtained with the SALT on 2013 April 27. 


which implies that MN18 is of spectral type of B1 or later 
(Walborn & Fitzpatrick 1990). 

Using the spectral atlas of OB stars by Walborn & Fitz¬ 
patrick (1990), we classify MN18 as a B1 la star. The la lu¬ 
minosity class also follows from the measured EW(H 7 ) of 
1.61±0.04 A and the calibrations of Balona & Crampton 
(1974; see their table V). 

The high interstellar extinction towards MN18 (see Sec- 
tions l4.2l and l5.l|) is manifested in numerous DIBs, the most 
prominent of which are at 4429, 4762, 4881, 5780, 5797, 6203, 
6269, 6283 and 6614 A. 

Inspection of the spectra of MN18 obtained in 2009- 
2014 did not reveal significant changes in their appearance, 
except for the variability of the Ha line. It was in emission 
in 2009 and 2013, in absorption in 2012, and in transition 
between these two states in 2014, with the line centre in ab¬ 
sorption and the line wings in emission. Comparison of the 
SALT high-resolution spectra obtained during two consec¬ 
utive nights in 2013 shows that the Ha line (and the C ll 
AA6578, 6583 doublet as well) changes its shape on a daily 
time-scale (see Fig. 51 ), which is typical of early-type super¬ 
giants (e.g. Morel et al. 2004; Martins et al. 2015). This vari¬ 
ability introduces an ambiguity in mass-loss estimates based 
on synthesizing the Ha line with model stellar atmosphere 
codes (see Section ing . 

Using the 1.9-m telescope and SALT spectra, we 
searched for shifts in RV as a result of possible binarity 
of MN18 and found that all measurements agree with each 
other within their error margins. 



Figure 4. Daily time-scale line profile variability in Ha and the 
Cll AA6578, 6583 doublet: 2013 April 26 (solid line) and 2013 
April 27 (dashed line). 












Wavelength (A) 

Figure 5. Parts of the normalized FEROS (upper panel) and high resolution SALT (bottom panel) spectra of MN18 (shown by blue 
solid lines), compared with the best-fitting model spectrum (red dotted line) with the parameters as given in Table[5] 



log Ilk 

Figure 6 . Observed flux distribution of MN18 in absolute units, based on the photometric measurements (blue rectangles) compiled in 
Tablcfn compared to the emergent flux of the reddened model spectrum (red dashed line) with the parameters as given in Tablc[3 


4.2 Spectral modelling 

For the spectral analysis we used spectra from three epochs: 
the low- and high-resolution SALT spectra taken on 2012 
March 2 and 2013 April 27, respectively, and the echelle 
FEROS spectrum taken on 2014 March 9. Our analysis is 
based on synthetic spectra computed with the non-LTE stel¬ 
lar atmosphere code CMFGEN (Hillier & Miller 1998). The 
following model atoms and ions were taken into account: H l, 
Hei-ii, Ci-iv, Niii-v, Oi-v, Mgii, Alii, Siii-iv and Fei-v. 
Because we are mainly interested in the evolutionary status 
of MN18, we varied only the abundances of He, C, N and O, 
while the abundances of other metals are set at solar values 


(taken from Asplund et al. 2005). Since our spectra did not 
cover the UV resonance lines, we cannot resolve the degen¬ 
eracy between the wind velocity law, /3, and a wind volume 
filling factor of optically thin (micro) clumps, /v. Similar 
line profiles for Ha can be achieved by varying and/or 
/v. Based on our model computations, j3 ranges between 1.5 
and 2.5, and fv between 0.25 and 1.0. In order to estimate 
the error bars in the effective temperature, Teg, mass-loss 
rate, M, and surface gravity, logg, we computed a grid of 
models around the preferred values. The best-fitting model 
is shown in Fig. [5] 

Teff is based on the ratio of the Fe ll and ill lines and the 
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intensity of the Hel, Nil and Sill ones, and is 21.1±1.0kK. 
We employed H/3 to determine logg=2.75±0.20. The uncer¬ 
tainty of log g is rather large as a result of low signal to noise 
ratio in the FEROS and low resolution in the SALT spectra. 
The terminal velocity, Voo , is difficult to constrain dne to the 
weak emission of Ha. No other suitable lines are available in 
the covered wavelength range. Therefore, the uncertainty is 
large and we estimate Vcx of Ril400±500 kms“^. For com¬ 
parison, we used the relation between Voo and the photo- 
spheric escape velocity, Wesc, from Tamers, Snow & Lind- 
holm (1995), which for Tefr = 21.1 kK, loggi=2.75 and the 
bolometric luminosity of log(L/ Lq)=5.42 (see below) gives 
Voo/vesc = 2.6 and Uck) = 1370±250 kms“^. 

The mass-loss rates for the three epochs were derived 
from Ha. Over these epochs the profile of Ha changed and, 
correspondingly, M as well. The obtained values for the 
individual epochs with log(L/ Lq)=5.42, Uoo = 1400 kms“], 
/3=2.0 and /v = 1.0 are log(M/ Mq yr-i)=-6.38±0.2^ 
(2012 March 2), -6.17T0.25 (2013 April 27), and 

—6.37T0.25 (2014 March 9). We observe a variation of M by 
a factor of Ril. 6 . The errors in M account for the uncertain¬ 
ties in log(L/L 0 ), Voo, /3-parameter and micro-clumping. 
Petrov et al. (2014) found that the optical depth of the 
Ha line below the bi-stability jump (Ten r;22.5 kK) increases 
with decreasing Teft. Correspondingly, the profile of this line 
becomes sensitive to macro-clumping and the derived M 
could be underestimated. However, the effective tempera¬ 
ture of MN18 is close to that of the bi-stability jump, while 
the Voo/veso ratio of 2.6 suggests that MN18 has wind prop¬ 
erties expected for a star on the hot side of this jump. 
Therefore, it appears likely that the influence of the macro- 
clumping on the mass loss rate of MN18 is small, i.e., negli¬ 
gible in comparison to the overall uncertainties. 

The CNO abundances are based on the ratio of the C, 
N and O spectral lines. Tabled shows that C is reduced 
by a factor of «3.5 and N is enriched by «4 times due to 
the CNO cycle. The abundance of O is still close to its ini¬ 
tial value within the error margins. At the given wavelength 
range the H and He lines are insensitive to changes in the 
He abundance, so that the He mass fraction could not be 
derived. 

To estimate the reddening towards MN18, we match 
the model spectral energy distribution (SED) with the ob¬ 
served one in the B, V and Ka bands (see Table 1). We ex¬ 
tracted intrinsic B — V and V — Ks colours from the model 
flux by applying approximate filter function^ and match 
the resulting values of the colour excesses E{B — V) and 
E{V — Ka) as described by Bestenlehner et al. (2011) to 
determine the total-to-selective absorption ratio Ry and V 
band extinction Ay. We derive log(L/ L0)=5.42±O.2 adopt¬ 
ing the absolute visual magnitude of Mv =—6.8T0.5 typical 
of Bla supergiants (Crowther et al. 2006). Accounting for 
the uncertainty in Tefr and photometric errors, the uncer¬ 
tainty in log(L/L 0 ) sums up to ±0.3 dex. Even though we 
find a variability in Ha, the temperature and bolometric cor- 


® This is an upper limit because the Ha line was in absorption 
in 2012. 

® The filter functions were approximated for each filter system in¬ 
dividually using Gaussian and boxcar functions with appropriate 
widths and wavelengths. 


Table 5. Stellar and extinction parameters for MN18. 


Teff (kK) 

21 . 1 T 1.0 

log(L/L 0 ) 

5.42T0.30 

R* (R©) 

on C- + 18.5 

oo.a_j^2,7 

logs 

2.75T0.20 

log(M/M 0 yr-l)(“) 

-6.38T0.25 

log(M/M 0 yr-l) 6 ") 

-6.17T0.25 

log(M/M 0 yr-l)(“) 

-6.37T0.25 

0 

1.5 to 2.5 

fv 

0.25 to 1.0 

Voo (kms“l) 

1400T500 

Ueq (kms“l) 

104±7 

C (mass fraction) 

6.5±2.5xl0-^ 

N (mass fraction) 

2.4j:^-3 X 10-3 

O (mass fraction) 

5.4^2 X 10-3 

E{B - U) (mag) 

1.97T0.10 

Rv 

3.26T0.10 

Mv (mag) (adopted) 

-6.8T0.5 

d (kpc) 

5.6tl;3 


Note: “2012 March 2; '’2013 April 27; “2014 March 9. 


rection (BC=—2.01 ma£l) are nearly unchanged, so that we 
can assume no significant variations in the luminosity over 
time. The resulting SED well matches the photometry from 
optical to mid-infrared wavelengths, and no infrared excess 
is observed (see Fig.|5]). 

By converting the luminosity range into distances, we 
found that MN18 lies at a distance of l kpc. The 

line of sight towards MN18 (ZRi321°) first crosses the Carina- 
Sagittarius Arm at «0.7-1.5 kpc and then the Crux-Scutum 
one at «3.6-4.4 kpc (Vallee 2014). Our distance estimate 
implies that MN18 is located either in the Crux-Scutum Arm 
or in the interarm region. In Section lO^ we argue that the 
former possibility is more plausible. Note that the luminosity 
of MN18 and its R* and M (given in Table[5]) can be scaled 
to different distances as <xd^, ocd and ocd^^^, respectively 
(Schmutz, Hamann & Wessolowski 1989). 

Using the FEROS spectra, we estimated the pro¬ 
jected rotational velocity, vsini, of MN18. For this, 
we used the relations between the FWHMs of the Hel 
AA4387, 4471 lines and wsini given in Steele, Negueru- 
ela & Clark (1999). For FWHM(4388)=2.24±0.25 A 
and FWHM(4471)=2.10±0.13 A, we found 
usinzRi94.1±10.5 kms“^ and 86.6T5.4 kms“^, respec¬ 
tively, with a mean value of 90.4T5.9 kms“^, which for 
i=60° (see Section[21) translates into the equatorial ro¬ 
tational velocity of iieq=104±7 kms“^. These estimates 
should be considered as an upper limit to the rotational 
velocity because the macroturbulence could contribute to 
the line broadening. 

A summary of the stellar and extinction parameters is 
given in TableO 


^ Note that almost the same value of BC (—2.03 mag) follows 
from relation between BC and Tgff of B supergiants given in 
Crowther et al. (2006). 
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Figure 7. A part of the 2D high-resolution spectrum of the equa¬ 
torial ring around MN18, showing nebular emission lines (from 
left to right) of [Nil] A6548, Ha, [Nil] A6548 and [Sii] AA6717, 
6731. 


5 SPECTROSCOPY OF THE EQUATORIAL 

RING 

5.1 Physical conditions and elemental abundances 

As noted in Section \n\ all SALT spectra were obtained 
with the slit oriented along the major axis of the equatorial 
nebula. In two-dimensional (2D) spectra the nebula appears 
as emission lines of H/3, Ha, [N ii] AA6548, 6584 and [S ii] 
AA6717, 6731, visible on both sides of MN18. Note that in 
the region of the Ha line the nebula is more than 5 mag 
fainter than the star and therefore it does not affect the 
stellar photometry appreciably. A part of one of the two 
2D high-resolution spectra is shown in Fig.[3 Besides the 
nebular spectrum and the spectrum of MN18, Fig. [3 also 
shows spectra of two other stars on the RSS slit, which are 
located at almost equal distances on both sides of MN18. In 
Section l^^ we show that the brighter star is of late 0-/early 
B-type. 

The upper and bottom panels of Fig.[8]plot the nebular 
emission line intensities and their ratios along the slit. From 
the upper panel it follows that the nebula can be traced up 
to ±11 arcsec (at a ten per cent level of the Ha peak inten¬ 
sity) from MN18, which corresponds with the extent of the 
equatorial ring in the Spitzer and SHS images (see Fig. [3). 
This panel also shows that the line intensities peak at «8.5 
arcsec from the star, which reflects the toroidal topology of 
the nebula. 

The lower panel of Fig.[8]shows that the [N ll] A6584/Ha 
line ratio increases with distance from MN18 and peaks 
at the edge of the nebula. Using the prescription by Do- 
pita (1973), one finds that [N ll]/Ha is proportional to 
exp( —21 855/Te), where Ue and Te are the electron 
number density and temperature, respectively. Thus, the be¬ 
haviour of [N ii] A6584/Ha could be explained (at least in 
part) by the temperature increase of the photoionized gas 
towards the outer edge of the nebula (caused by the harden¬ 
ing of the ionizing radiation field), provided that the density 
is constant throughout the nebula. 

One-dimensional spectra (ID) of the nebula were ex¬ 
tracted by summing up, without any weighting, all rows 
from the area of an annulus with an outer radius of 10 arc¬ 
sec centred on MN18 and the central ±3 arcsec excluded. 
Both low- and high-resolution ID spectra of the nebula are 


Table 6. Line intensities of the MN 18 nebula. 


Ao(A) Ion 

F(A)/F(H/1) 

I(A)/I(H/3) 

4861 np 

1.000T0.236 

1.000T0.249 

6548 [N ii] 

7.083T1.183 

0.578T0.111 

6563 Ha 

36.750T6.136 

2.947T0.564 

6584 [N ii] 

21.792T3.647 

1.705T0.327 

6717 [S ii] 

3.708T0.631 

0.248T0.049 

6731 [S ii] 

3.750T0.637 

0.247T0.049 

C(H/3) 

3.30T0.22 

E{B-V) 

2.24T0.15 mag 

ne([Sll ]) 

560t«o 

_Q 

cm 


presented in Fig. [9] The emission lines detected in the spec¬ 
tra were measured using the programs described in Kni¬ 
azev et al. (2004). Table|6] lists the observed intensities of 
these lines normalized to H/3, F(A)/F(H/3), the reddening- 
corrected line intensity ratios, I(A)/I(H(1), and the loga¬ 
rithmic extinction coefficient, C'(H(1), which corresponds to 
E{B — U)=2.24±0.15 mag. In Table[6]we also give the elec¬ 
tron number density derived from the intensity ratio of the 
[Sii] AA6716, 6731 lines, ne([Sii]). The obtained value is 
comparable to densities derived for equatorial rings of other 
bipolar nebulae around blue supergiants and cLBVs (see 
Section[6]2]) . Both (^(H/J) and ne([Sll]) were calculated un¬ 
der the assumption that Te = 10"^ K. For lower temper¬ 
atures their values do not change much, e.g. C'(H/3)=3.27 
and rie = 510 cm“® if Te=7000 K. 

The electron number density could also be derived from 
the surface brightness of the equatorial ring in the Ha line, 
Sna, measured on the SHS image (cf. equation (4) in Frew 
et al. 2014): 



where I is the line of sight thickness of the ring and 
1 R=1 Rayleigh=5.66x 10~^® erg cm“^ s~^ arcsec”^ at Ha 
(here we assume that rie is constant within the ring and 
Te = 10'* K). 

Using equations (1) and (2) in Frew et al. (2014) and the 
flux calibration factor of 7.2 counts pixeU^ R“^ from their 
table 1, and adopting the observed [N ll] to Ha line intensity 
ratio of 0.79 from Table[5](here [N ll] corresponds to the sum 
of the A6548 and A6584 lines), we obtained the peak surface 
brightness of the northeastern (brightest) half of the ring 
(corrected for the contribution from the contaminant [N ll] 
lines) of 5 'hq~ 120 R. Using equation (1) with E{B — V)=1.97 
mag and Z«0.47 pc (the line of sight thickness of the torus), 
one finds ne~570 cm~®, which agrees well with ne([Sll]) 
given in Table[6l 

To obtain an order of magnitude estimate of the mass 
of the equatorial ring, ATing, we assume that the ring has a 
toroidal geometry with the major and minor radii of R—6.5 
arcsec (r; 0.17 pc) and r=4.5 arcsec («0.12 pc), respectively, 
and that its material is fully ionized and homogeneous. In 
doing so, one finds Mring~l Mq. The actual value of Mring 
could be higher if the ionized gas encloses neutral and, cor¬ 
respondingly, more dense material. 

The [N ii] and [S ll] line intensities can be used to 
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Figure 8. Variation of emission line relative fluxes and their ra¬ 
tios along the slit. The solid vertical line corresponds to the po¬ 
sition of MN18, while vertical dashed lines at ±2.5'' from the 
solid line mark the area where no nebular emission was detected. 
NE—SW direction of the slit is shown. Upper panel: The rela¬ 
tive fluxes of the Ha, [Nil] A6584 and [Sii] A6717±A6731 lines 
after continuum subtraction. Bottom panel: [Nil] A6584/Ha and 
[Sii] A6717±A6731/[N ii] A6584 line ratios. 



Figure 9. ID low-resolution spectrum of the equatorial nebula 
around MN18. The insert shows a part of the high-resolution 
spectrum of the nebula taken on 2013 April 27. 


estimate the nitrogen to sulphur abundance ratio, which 
is almost independent of rie and Te, provided that rie ^ 
1000 cm~^ and Te ^ 10^ K. In this case, the abundance ra¬ 
tio is given by (see Benvenuti, D’Odorico & Peimbert 1973 
and references therein): 


A^(N+) ^ 7(6584) 

iV(S+) 7(6716 + 6731)' 

Using equation ([2]) and Tablejd] 


( 2 ) 


one finds 



Radius (arcsec) 


Figure 10. Upper panel: Ho heliocentric radial velocity distribu¬ 
tion along the slit. The solid vertical line corresponds to the po¬ 
sition of MN18, while the vertical dashed lines at ±2.5" from the 
solid line mark the area where no nebular emission was detected. 
The model fit is shown as a (blue) solid line. The (red) dashed and 
dotted horizontal lines show, respectively, the systemic velocity 
of the equatorial ring and its error margins. Bottom panel: The 
difference between the observed and model radial velocity distri¬ 
butions. The mean residual velocity and its errors are shown with 
the (blue) solid and dotted lines, respectively. 


Ar(N+)/7V(S+)«12.43l®;°i, which is «2.4lJ;2 times 
larger than the solar value of 5.12 (Asplund et al. 2009). 
Since the sulfur abundance is not expected to change 
during stellar evolution, one can argue that the high N to 
S ratio exists because of an elevated nitrogen abundance 
in the line-emitting material. Thus, using the solar N and 
S abundances from Asplund et al. (2009), one finds a N 
abundance of the ring of log(N/H) + 12=8.21. 


5.2 The position—velocity diagram and kinematics 
of the equatorial ring 

Fig.[T0] plots the distribution of the heliocentric radial ve¬ 
locity of the Ha line, Uhei, along the slit, calculated using 
the method and programs described in Zasov et al. (2000). 
The systematic shifts originating from the RSS flexure were 
excluded by taking into account the closest bright night sky 
line. Finally, only those velocity measurements were used 
which satisfy the criteria S/N > 3 and (t„ < 6 kms“^. 

To estimate the systemic velocity of the equatorial ring, 
Vaya, and its expansion velocity, Uexp, we fit the observed 
distribution of Uhei with the model: 

Uhel — Uexp COS a + Vaya , (3) 

where Uexp and Usys are free parameters, and a is the an¬ 
gle between the line of sight and Uexp at each point of the 
ring. The result of the fit is shown as a (blue) solid line 
in the upper panel of Fig. HOI The (red) dashed and dot¬ 
ted horizontal lines show, respectively, the derived Ways of 
—66.0 kms“^ and its error margins (±1.3 kms“^). We also 
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found Vexp = 7.9±1.2 kms“^. The difference between the 
observed and model velocity distributions (the residual ve¬ 
locity hereafter) is shown in the bottom panel of Fig. llOi 
where the (blue) solid and dotted lines indicate, respectively, 
the mean residual velocity of —0.05 kms“^ and ±1(7 of the 
mean (±1.8 kms“^). The residual velocity distribution does 
not show any systematics, which allows us to set an upper 
limit of 1.8 kms“^ on the rotational velocity of the ring. 

The radius of the ring of r—0.29 pc together with Vexp 
yield the kinematic age of the ring of tkin=r/uexpSs3.7 x 10^ 
yr. Assuming that the equatorial ring and bipolar lobes are 
the result of the same episode of enhanced mass loss from 
MN18 (i.e. their tkin are the same; cf. Brandner et al. 1996b), 
one finds that the lobes should expand at «60 kms“^. Us¬ 
ing Mring from Section lSm and tkin, one can also find an 
order of magnitude “kinematic” estimate of M during this 
episode, i.e. Mkin'~iV/'ring/tkin«2.7xlO“® MQyr“^, which is 
about two orders of magnitude higher than the present-day 
values (see Table[5l). The discrepancy between the mass-loss 
rates would be somewhat smaller if the expansion of the 
equatorial ring was faster initially. (Correspondingly, the ex¬ 
pansion velocity of the lobes would be smaller as well.) On 
the other hand, the duration of the high mass loss episode 
could be shorter than fkin, which would result in an even 
higher Mkin- 


6 DISCUSSION 

In Section[Tl we noted the existence of a class of Galactic and 
extragalactic hourglass-like nebulae associated with blue su¬ 
pergiants and mentioned seven known members of this class. 
One of them -- the three-ring nebula associated with the blue 
supergiant Sk—69°202 (the progenitor star of the SN 1987A) 
- was formed shortly (~10^ yr) before the star exploded as 
a supernova. Currently, one can observe how the supernova 
blast wave interacts with the equatorial ring of the nebula 
(Sonneborn et al. 1998; Sugerman et al 2002; Groningsson et 
al. 2008; Fransson et al. 2015). It is also believed that LBVs 
could be immediate precursors of supernova explosions (e.g. 
Kotak & Vink 2006; Groh et al. 2013; Justham et al. 2014). 
Thus, one can expect that some other massive stars with 
bipolar nebulae will soon explode as supernovae and that 
the appearance of their young (meaning from decades to 
centuries old) diffuse supernova remnants would be affected 
by the dense material of the nebulae. It is therefore worth 
to compare the derived parameters of MN18 with those pre¬ 
dicted by stellar evolutionary models (e.g., Brott et al. 2011; 
Ekstrom et al. 2012) to determine the evolutionary status 
of this star and thereby to understand whether or not it is 
close to explode as well. Knowing the evolutionary status of 
MN18 would also allow us to constrain the mechanism for 
the origin of the bipolar nebula around this star, which in 
turn would help to better understand the origin of bipolar 
nebulae around other massive stars. 

6.1 Evolutionary status of MN18 

The range of luminosities derived for MN18 in Section lT^ 
corresponds to the initial mass, Minit, of this star of 
«30±10 Mq. Using the grid of rotating main sequence stars 
by Brott et al. (2011), we found that the main parameters 


of MN18 (temperature, luminosity, gravity and CNO abun¬ 
dances) can be matched very well with models with Minit 
in the above range, provided that the initial rotational ve¬ 
locity of the star, uinit, was of ~200-400 kms“^. For ex¬ 
ample, the model of a 30 Mq star with Uinit = 321 kms“^ 
predicts that after «5.5 Myr the star would have the 
following parameters: Teff^Sl.OTO.S kK, log(Z//L0)«5.41, 
i?,«38±2 Rq, logg«2.70±0.05 and a current rotational ve¬ 
locity of «110±10 kms“^. The CNO abundances of this 
star would be «6.3 x 10“"', 1.6 x 10“^ and 3.4 x 10“®, 
respectively. All these parameters are in good agreement 
with those derived for MN18 (see Table[3. The models with 
low Uinit also match very well the Tefi, log(Z//LQ) and logg 
values of MN18, but their CNO abundances differ signif¬ 
icantly. Specifically, the nitrogen abundances predicted by 
these models are much lower than what follows from our 
spectral modelling. 

Similarly, we found that the rotating models by Ek¬ 
strom et al. (2012) also fit well the parameters of MN18 
provided that this star just left the main sequence, while 
in the case of a star undergoing a blue loop evolution the 
CNO abundances predicted by the models differ significantly 
from those of MN18. The non-rotating models by Ekstrom 
et al. (2012) do not fit the CNO abundances of MN18 as 
well. Proceeding from this, one might conclude that MN18 
is a redward evolving star, which has only recently become 
a blue supergiant (cf. Gvaramadze et al. 2014a). This con¬ 
clusion agrees with the suggestion by Earners et al. (2001) 
that LBV nebulae could be formed already at the end of 
the main sequence phase and that their enhanced nitrogen 
abundances are due to rotationally induced mixing in their 
underlying stars. Similar suggestions were put forward by 
Smartt et al. (2002) and Hendry et al. (2008) to explain the 
CNO abundances of the bipolar nebula around Sher25. 

The duplicity of massive stars provides another clue for 
understanding the origin of enhanced N abundance in main 
sequence stars because accretion of CNO-processed material 
from the mass donor star or merger of the binary compo¬ 
nents can enhance the surface N abundance of the mass 
gainer star or the merger product, respectively, by a factor 
of several even in the absence of rotational mixingj (Langer 
2012; Glebbeek et al., 2013). As follows from the studies 
by Mason et al. (2009), Sana & Evans (2011) and Chini 
et al. (2012), the majority of massive stars are members of 
close binary systems. Furthermore, according to Sana et al. 
(2012), about 70 per cent of all massive stars are strongly 
affected by binary interaction during their lifetimes through 
the mass transfer, common envelope evolution or merger, 
and for more than half of them the interaction takes place 
already during the main sequence phase (see also de Mink et 
al. 2014). All of these types of binary interaction can enhance 
the surface N abundance already during the main sequence 
phase, and are also the likely culprits of the production of 
substantial mass loss in the equatorial plane. 

Thus, MN18 could be a mass gainer from a Case A or 
Case B system, and as such it could either be rejuvenated 
or not. In the first case there would be little difference to 


® Models taking into account rotational mixing in the mass gain¬ 
ers show that the N abundance can additionally be enhanced by 
a factor of 3 (Langer et al. 2008). 
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single star evolution, while In the second one the star at 
Its position in the Hertzsprung-Russell diagram could either 
be core hydrogen or core helium burning (Braun & Langer 
1995; Justham et al. 2014). The latter possibility, however, 
is less likely because the N abundance of the star would be 
much higher than what we found for MN18. 

To summarize, the above considerations suggest that 
MN18 is near the end of the main sequence phase, which 
makes it unlikely that this star will soon explode as a super¬ 
nova. 

6.2 Origin of the bipolar nebulae 

Now we discuss the origin of bipolar nebulae around MN18 
and other massive stars. The rotational symmetry of these 
nebulae implies that the stellar rotation plays a significant, 
if not a decisive, role in their origin. 

6.2.1 Single star scenarios 

Tamers et al. (2001) proposed that at the end of the main 
sequence the initially rapidly rotating stars may reach the 
limit of hydrostatic stability (fl-limit; Langer 1997, 1998) 
in their outer equatorial layers, where the joint action of 
the centrifugal force and the radiation pressure balances the 
gravity (cf. Langer, Garia-Segura & Mac Low 1999). This 
may result in ejection of a large amount of mass (up to 
several M©) mainly in the equatorial plane. The ejection 
process could be recurrent with time intervals of ~ 10® —10^ 
yr, which are related to the thermal time-scale of the stars. 
Correspondingly, the stellar wind in between the ejections 
could be collimated in the polar direction by the presence of 
the dense equatorial material. 

To reach the fl-limit the star should be luminous enough 
to arrive at or exceed its Eddington limit in the course of 
its lifetime and/or be a rapid rotator. For single stars the 
first condition could be satisfied during the advanced evolu¬ 
tionary phases of very massive (>50 Mq) stars when their 
masses are considerably reduced due to the extensive stellar 
wind. Such stars can achieve the fl-limit even if they are 
slow rotators (e.g. Langer & Heger 1998). The resulting en¬ 
hanced mass loss of these stars is almost spherically symmet¬ 
ric and, correspondingly, the circumstellar nebulae around 
them should be almost spherical as well. In the case of very 
rapid rotators, the fl-limit could be reached already during 
the main sequence. This is conceivable even for those stars 
which never approach the Eddington limit during their life¬ 
time. Correspondingly, the mass loss from these stars would 
be highly aspherical, with a slow and dense wind at lati¬ 
tudes close to the equator and a fast and more tenuous one 
in the polar direction. These stars would produce elongated 
or bipolar nebulae. 

Interestingly, the Brott et al.’s models predict that a 
20—30 Mq star with Uinit^SOO—600 kms“^ can achieve the 
fl-limit by the end of the main sequence, which results in a 
strong mass loss with M of >10“^ Mq yr“^. But the higher 
Uinit, the stronger the mixing, the longer the duration of 
the main sequence phase and the bluer the star. Corre¬ 
spondingly, it takes longer for the star’s Tes to decrease 
to 21 kK. The prolonged rotational mixing in turn results 
in CNO abundances much different from those derived from 
our spectral analysis. 


The question on whether the surface velocity of single 
stars may hit critical rotation during the main sequence evo¬ 
lution was studied in detail by Ekstrom et al. (2008). They 
computed stellar models of various Minit, Uinit and metallic- 
ities, and found that the fl-limit could be reached by B stars 
with Uinit in the upper tail of the initial rotational velocity 
distribution. This could be responsible for the origin of the 
Be phenomenon. Ekstrom et al. (2008) also found that much 
greater mass loss from more massive stars at high (e.g. solar) 
metallicities prevents these stars from reaching the critical 
rotation. 

Heger & Langer (1998) proposed a mechanism which 
may spin-up the surface layers of rotating single stars to 
critical rotation after they leave the red supergiant branch. 
It was suggested that this mechanism may result in signifi¬ 
cant rotation of the progenitor star of the SN1987A (which 
performed a blue loop just before the supernova explosion) 
and thereby in the origin of the equatorial ring around it. 
The CNO abundances of the material ejected by a post- 
red supergiant star, however, would significantly differ from 
those of MN18 (see Section ISTII) . 

Another possibility to increase mass loss rates of sin¬ 
gle stars is based on the bi-stability in the radiation driven 
winds of supergiants near the spectral type B1 (Pauldrach & 
Puls 1990; Tamers & Pauldrach 1991), which is manifested 
in a drastic change in the ionization states of Fe when Tefi 
drops a little below some critical value («21-25 kK; Vink, 
de Koter & Tamers 1999). This change leads to a factor 
of 2 decrease in Voo, from iioo/uesc~2.6 for supergiants of 
types earlier than B1 to Uoo/uesc~1.3 for B supergiants of 
later types (Tamers et al. 1995), and in a factor of several 
increase in M (Tamers et al. 1995; Vink et al. 1999). The 
increase of M caused by the bi-stability jump could be re¬ 
sponsible for the ejection of shells by the prototype LBV 
P Cygni (Pauldrach & Puls 1990) and the observed steep 
drop in the rotational velocities of massive stars at Tefj of 
22 kK (Vink et al. 2010). 

More importantly, the enhanced mass loss caused by 
the bi-stability mechanism could be asymmetric if the star 
is a rapid rotator. In this case, the gravity darkening (the 
von Zeipel effect) leads to a lower Tefr at the equator (e.g. 
Maeder 1999) and the bi-stability jump first occurs near 
the equator as well. This in turn results in enhanced mass 
loss at low stellar latitudes. Tamers & Pauldrach (1991) and 
Tamers et al. (1995) suggested that the bi-stability jump in 
the winds of rapidly rotating early B-type supergiants could 
be responsible for the origin of dense equatorial outflows 
(excretion disks) around B[e] stars. The detailed study of 
this possibility by Pelupessy, Tamers & Vink (2000), how¬ 
ever, showed that the bi-stability mechanism can produce 
only about a four-fold increase in M, which is not enough 
to explain the equator to pole density contrast of a factor 
of «100 observed for the winds of B[e] stars. It is therefore 
doubtful that the bi-stability jump alone could explain the 
about two order of magnitude difference between M derived 
from our spectral modelling of MN18 (see Table[5]) and M 
required to produce the equatorial ring around this star (see 
Section lT^ . 

To conclude, although we cannot exclude the possibility 
that some single massive stars can hit critical rotation in 
the course of their evolution (resulting in massive equatorial 
outflows), we believe that the origin of the equatorial ring 


around MN18 is more likely because this star is or was a 
member of a close binary system. 


6.2.2 Binary star scenarios 

The binary population synthesis modelling by de Mink et 
al. (2013) suggests that about 20 per cent of massive main 
sequence stars should be rapid (> 200 kms“^) rotators be¬ 
cause of mass transfer and mergers in close binary systems. 
About 2 per cent of these stars are expected to have very 
high (> 400 kms“^) rotational velocities (with the most 
extreme cases reaching «600 kms“^) of which a quarter 
are the merger remnants, de Mink et al. (2013) noted that 
their predicted distribution of rotational velocities of main 
sequence stars agrees well with the observed one, which 
points to the possibility that mass transfer and mergers are 
the main causes for the rapid rotation of massive stars (cf. 
Packet 1981; Langer et al. 2003; Dufton et al. 2011). It is 
tempting therefore to speculate that these two binary inter¬ 
action processes are the main causes for the origin of nearly 
critically rotating stars and equatorial outflows from them. 

The mass transfer in binary systems may not only spin- 
up the mass gainers, but under proper conditions it could 
also result in the origin of circumbinary disk-like structures 
by means of the common envelope ejection (e.g. Morris 
1981), mass loss through the second Lagrange point (e.g. 
Livio, Salzman & Shaviv 1979; Petrovic, Langer & van der 
Hucht 2005; Smith et al. 2011), and focusing of the primary’s 
wind towards the orbital plane by the gravitational field of 
the companion star (e.g. Fabian & Hansen 1979; Mastrode- 
mos & Morris 1999). Similarly, the binary mergers can not 
only produce the rapidly rotating stars, but also result in 
mass ejection (~10 per cent of the binary mass; Suzuki et al. 
2007), which is mainly concentrated in the binary’s equato¬ 
rial plane. Moreover, the instant mass loss during the merger 
event may bring the merger product to the Eddington limit, 
which in turn would result in a copious mass loss by stellar 
wind. 

If the origin of the equatorial ring around MN18 is in¬ 
deed caused by the equatorial mass loss from either a nearly 
critically rotating mass gainer or a merger remnant, then the 
present moderate rotational velocity of this star would im¬ 
ply that the star has suffered significant angular momentum 
loss on a time-scale of tkin- The most efficient mechanisms 
for the spin down of (single) massive stars is the magnetic 
braking via a magnetized wind (Poe, Friend & Cassinelli 
1989; ud-Doula, Owocki & Townsend 2009), which for stel¬ 
lar magnetic fields of ~10® — 10^ G and sufficiently high M 
can result in spin-down time-scales as short as ~10"^yr. It 
is not clear whether initially single stars can simultaneously 
possess such high magnetic fields and be rapid rotators, but 
it is believed that very strong large-scale magnetic fields 
can intermittently be created because of the strong shear 
in merging binary systems (Langer 2012; Wickramasinghe, 
Tout & Ferrario 2014). Such magnetic fields can brake the 
rapid rotation of the merger remnants through dynamical 
mass loss. Polarimetric observations of MN18 and other stars 
with bipolar nebulae would be of great importance to check 
whether they possess large-scale magnetic fields of strength 
of ~100-1000 G, which are not typical of the majority of 
massive main sequence stars (e.g. Wade et al. 2014; Morel 
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et al. 2015) and whose origin can be the consequence of bi¬ 
nary merger. 

In Section[S31 we assumed that the kinematic ages of 
the equatorial ring and bipolar lobes of the nebula around 
MN18 are the same (i.e. both components of the nebula 
were formed simultaneously). But it is quite possible that 
the ring was formed in the first instance (e.g. because of 
a dynamical mass ejection caused by the binary merger) 
and that later on its dense material collimated the stel¬ 
lar wind into bipolar lobes. If so, then pure ring-like neb¬ 
ulae should exist whose central stars still have had no 
time to blow-up pronounced bipolar lobes and which still 
might be fast rotators. Indeed, several examples of such 
rings exist. One of them, named [GKF2010] MN109 in the 
SIMBAD data base, encircles the candidate blue super¬ 
giant 2MASS J19281457+1716230 (Phillips & Ramos-Larios 
2008). This ring is best visible at 8/rm (see the right-hand 
panel of Fig. [13, while at 24 /rm there are no signatures of 
an extended bipolar outflow. Another example, known as 
[GKF2010] MN56, is associated with the cLBV MN56. In 
Fig. ll2l we present for the first time the AT-band spectrum 
of this star obtained with the ISAAC spectrograph on the 
Very Large Telescope (VLT) during our observing run in 
201(0. The spectrum shows strong emission lines of H, Hel, 
Fell and Mgll and is almost identical to that of the cLBV 
HD 316285 (see fig. 2 in Hillier et al. 1998). [GKF2010] MN56 
also has a clear-cut ring-like morphology with a “horn” at¬ 
tached to the south-east edge of the ring (see the right-hand 
panel of fig. 4 in Paper I for the IRAC 5.8 ^m image of this 
structure). The presence of the horn suggests that the stel¬ 
lar wind already starts to be collimated into bipolar lobes. 
The saturated 24 /rm image of [GKF2010] MN56 (see the 
left-hand panel of fig. 4 in Paper I) did not allow us to dis¬ 
cern the horn, but it also did not show signatures of a more 
extended bipolar outflow. High-resolution spectroscopy of 
2MASS J19281457-I-1716230 and MN56 would be desirable 
to determine the rotational velocities of these two stars, as 
well as the evolutionary status of the former one. 

Binary population studies suggest that the rate of mas¬ 
sive binary mergers in the Galaxy is ~10-20 per cent of the 
Galactic supernova rate, which corresponds to one per 200 
yr (Langer 2012). Assuming that all mergers result in the 
formation of ring-like or bipolar nebulae and that the disper¬ 
sion time-scale of these nebulae is ~10^ yr, one finds that our 
Galaxy should contain about 100 such nebulae. The actual 
number of these nebulae however should be much smaller 
because the majority («80 per cent; Gies 1987) of massive 
stars reside in their parent star clusters where the stellar 
winds and ionizing emission of neighbouring hot massive 
stars could prevent the formation of discernible circumstel- 
lar nebulae. This inference could be supported by the fact 
that all massive stars with known well-shaped (bipolar or 
ring-like) nebulae are located either at the periphery of their 
parent star clusters or in the field, and therefore are likely 
runaway stars (Gvaramadze et al. 2012a,b; 2014b). 


® Full details of these observations will be presented elsewhere, 
while some preliminary results could be found in Stringfellow et 
al. (2012a,b). 
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Figure 11. Spitzer MIPS 24/jm (left) and IRAC 8/^m images of a ring-like nebula around the possible blue supergiant 
2MASS J19281457+1716230. 



Figure 12. Normalized VLT/ISAAC A-band spectrum of the 
cLBV MN56 with principal lines indicated. 

6.3 MN18 and other B supergiants with bipolar 
nebulae 

All but one (MN13) of seven hourglass-like nebulae around B 
supergiants mentioned in Section[T] are visible (at least par¬ 
tially) at optical wavelengths, and five of them were studied 
by means of optical spectroscop\F^. In Table[7l we compare 
the main parameters of the central stars and equatorial rings 
of these five nebulae. 

TableQ shows that the luminosities of the stars as¬ 
sociated with bipolar nebulae span a large range, from 
log(I//L 0 )= 4.7 to 5.8 (with the luminosity of MN18 be¬ 
ing close to the median value of «5.3). From this one can 

There are no reported spectroscopic observations of the equa¬ 
torial belt (visible in optical light; Marston & McCollum 2008) of 
the bipolar nebula around the cLBV MWC 349A. 


conclude that the stellar luminosity (or Minit) is not a deci¬ 
sive factor in formation of bipolar nebulae, which is expected 
if the nearly critical rotation is the main reason for the ori¬ 
gin of these nebulae. In this regard, it is worthy to note that 
some PNe also possess the hourglass-like morphology, whose 
origin is most likely due to binary interactions (e.g. Corradi 
& Schwarz 1993; Sahai et al. 1999; De Marco 2009). 

One can also see that Tefj values span quite a large 
range as well, from 13 to 22 kK, which might be considered 
as an indication that the bi-stability mechanism may not 
contribute to the origin of at least some of the bipolar neb¬ 
ulae. This inference, however, could be inconclusive because 
Tefi of bona fide LB Vs can vary in the same range on time 
scales as short as several years (e.g. Stahl et al. 2001; Groh 
et al. 2009). One cannot therefore exclude the possibility 
that Teff of Sk—69°202 and HD 168625 were close to the bi¬ 
stability limit during the formation of nebulae around these 
stars. 

The axial symmetry of the bipolar nebulae suggests that 
the tilt angles of their equatorial rings are equal to those of 
the rotational axes of the central stars. This allows one to 
estimate the rotational velocities for those stars whose pro¬ 
jected rotational velocities are known. Tabled shows that 
none of these four stars are rapid rotatorf^. Thus, if the 
origin of the equatorial rings is due to the nearly critical 
rotation of their underlying stars, then the spin-down time- 
scales of these stars should be shorter than tkin of the equa¬ 
torial rings, all of which are ^lO'^ yr. It is reasonable to 
associate this drastic spin down with a brief episode of very 
high M (e.g. Heger & Danger 1998), which is much higher 
than the current values of M and perhaps even higher than 
the “kinematic” ones, Mkin ~ ALring/tkin- 

Table[7]also shows that the linear sizes of all equatorial 
rings are almost the same: 0.2—0.3 pc. The largest value of 
0.3 pc was obtained for the equatorial ring around MN18 
under the assumption that the distance to this object is 
5.6 kpc. At a shorter distance of ~4A kpc (allowed by the 
error margins of d; see also Section lOTdll the radius of the 


There is no information on the projected rotational velocity 
of the progenitor star of the SN1987A. 
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ring would be similar to those of other equatorial rings. The 
similarity in size of the rings produced by stars of different 
Minit and evolutionary status, however, is most probably 
a coincidence, because it is difficult to imagine any viable 
mechanism for constraining it. 

The number densities of all equatorial rings are also 
about the same 1000 cm“®) except of the ring around 
Sk—69°202, for which an order of magnitude higher value 
was derived from the fading of the [O ill] A5007 emission 
line (Plait et al. 1995). The masses of the equatorial rings 
range from «0.1 to 1 M© with the lowest mass belonging to 
the ring around Sk—69°202. This ring, however, might be 
more massive if we see only its inner-faced ionized part (e.g. 
Plait et al. 1995; cf. Fransson et al. 2015). 

Finally, Table[7] gives also the nitrogen abundances (by 
number) of the equatorial rings. In the Sk—69°202 ring the 
N abundance is a factor of 20 greater than the LMC in 
general (7.14 dex; Russell & Dopita 1992), which is consis¬ 
tent with the fact that the material of this ring was ejected 
during the final stage of the evolution of Sk—69°202. The 
Sher25 ring is also greatly enriched with nitrogen as com¬ 
pared to the background H ll region NGC 3603 and the Sun, 
respectively, by factors of «30 and 12 (Hendry et al. 2008). 
Although this overabundance could be interpreted as an in¬ 
dication that Sher 25 have evolved through the red super¬ 
giant stage (Brandner et al. 1997a,b), Hendry et al. (2008) 
argued that Minit of Sher 25 of ~50 Mq rather implies that 
this star did not experience this stage at all. Instead, Hendry 
et al. (2008) suggested that the enhanced N abundance is 
because of rotationally induced mixing while the star was 
on the main sequence. Mass accretion or a merger in a bi¬ 
nary system might be responsible for this enhancement as 
well (cf. Section lSTTIl . The abundances of the rings produced 
by HD 168625 and MN18 exceed the solar abundance by a 
factor of few, while that of the [SBW2007] 1 ring is subso¬ 
lar. It is therefore likely that the underlying stars of these 
three rings have not yet passed through the red supergiant 
stage and that all these rings were formed at the end of the 
main sequence stage or soon after it (cf. Tamers et al. 2001; 
Smith 2007; Smith et al. 2007; Gvaramadze et al. 2014a). 
The largely different N abundances of the equatorial rings 
imply that blue supergiants are capable to produce them in 
very different evolutionary stages, from the main sequence 
stage up to the pre-supernova stage. MN18 adds to the ev¬ 
idence that multiple ring nebulae could be formed rather 
early during the stellar evolution. 

6.4 MN18 and LyngaS 

Is Section[2l we mentioned that Lynga (1964a) listed MN18 
among the members of the newly identified open star clus¬ 
ter LyngaS. Lynga (1964b) noted that the cluster is strongly 
reddened, which did not allow him to estimate the interstel¬ 
lar extinction and the distance to the cluster. 

A literature search revealed three subsequent studies of 
LyngaS, based either on infrared or optical photometry. Gar- 
raro et al. (2006) constructed a colour-magnitude diagram 
of what they assumed to be a cluster using their own GGD 
BVI photometry and did not find any signihcant difference 
with respect to a comparison field. This led them to con¬ 
clude that “LyngaS is not a cluster, but an overdensity of a 
few bright stars probably located inside the Garina arm”. It 


turns out, however, that they confused the 1950 epoch co¬ 
ordinates with the 2000 epoch ones, which resulted in more 
than 0?5 offset between the group of stars they analyzed and 
the centre of the cluster given in Lynga (1964a,b). Moreover, 
the inspection of the finding chart of LyngaS given in fig. 13 
of Lynga (1964a) revealed that the geometric centre of the 
putative cluster is offset by about 5!5 from the cluster’s po¬ 
sition listed in Lynga (1964a,b; see tables V in these papers). 
This explains why Garraro et al. (2006) came to the conclu¬ 
sion that Lynga 3 is not a true cluster. 

More recently, LyngaS was studied by Kharchenko et 
al. (2013) and Dias et al. (2014), who used photometric and 
proper motion measurements from different catalogues to 
evaluate the membership probabilities for individual stars 
in the cluster region and to derive the distance to and age of 
the cluster. Although the coordinates of the cluster given in 
these two references are in a quite good agreement with each 
other (as well as with the position of the cluster’s geometric 
centre implied by the finding chart in Lynga 1964a), the ob¬ 
tained distance estimates, 4.2 kpc (Kharchenko et al. 2013) 
and 1.4 kpc (Dias et al. 2014), differ from each other drasti¬ 
cally. Moreover, the age estimates of logt=8.92 (Kharchenko 
et al. 2013) and 8.55 (Dias et al. 2014) are much older than 
the age of MN18 of logt«7 (see Section lOTTIl . Since this huge 
age difference cannot be explained by rejuvenation of MN18 
because of mass transfer or merger in a binary system, it 
is reasonable to assume that MN18 and the putative star 
cluster are projected by chance along the same line-of-sight. 
This raises the question where the parent cluster of MN18 
is? 

A possible answer to this question is that MN18 is a 
runaway star, which was ejected and arrived at its cur¬ 
rent position from a nearby massive star cluster (cf. Gvara¬ 
madze, Pflamm-Altenburg & Kroupa 2011; Gvaramadze et 
al. 2012a). The runaway status of this star is suggested 
by proper motion measurements given in the SPM4.0 (Gi¬ 
rard et al. 2011) and UGAC4 (Zacharias et al. 2013) cat¬ 
alogues (see Table[8]). Using these measurements, Usys de¬ 
rived in Section lA^ the Solar galactocentric distance i?o = 
8.0 kpc and the circular Galactic rotation velocity 0o = 
240 kms“^ (Reid et al. 2009), and the solar peculiar motion 
{Uq,Vq,Wq) = (11.1,12.2,7.3) kms“^ (Schonrich, Binney 
& Dehnen 2010), we calculated the peculiar transverse (in 
Galactic coordinates) velocity, m and Vh, the peculiar radial 
velocity, Vi, and the total space velocity, w*, of MN18 (see 
Table[Sl). To take into account uncertainties in the distance 
estimate, we adopted three values of the distance: ci=4.4, 5.6 
and 7.1 kpc (see Section lT^ . For the error calculation, only 
the errors of the proper motion and the systemic velocity 
measurements were considered. 

From Table[5]it follows that if the SPM4.0 and UCAG4 
proper motion measurements are correct then MN18 is a 
high-velocity runaway star, moving almost along the Galac¬ 
tic plane and the plane of sky. In this case, the parent cluster 
of MN18 would be at ~15°-25° from the current position of 
the star (i.e. at Z«336°-346°), provided that it was ejected at 
the very beginning of the cluster’s dynamical evolution (i.e. 
~5-10 Myr ago). It should be noted, however, that although 
the high space velocity derived for MN18 is not unusual for 
single runaway OB stars, it would be unrealistically high 
if this star have been ejected as a binary (cf. Section lR^ . 
Moreover, the ram pressure of the interstellar medium due 
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Table 7. Blue supergiants and cLBVs with hourglass-like circumstellar nebulae. 



Sk-69°202 

Slier 25 

HD 168625 

[SBW2007] 1 

MN18 

Spectral type 

B3l(i) 

B1.5Iab(2) 

B 6 Iap(®) 

BlIab(V 

B1 la 

log(L/L 0 ) 

^5(5) 

5.8(®) 

5.0-5.4(1) 

4.7(8) 

5.4 

Teff (kK) 

16(®) 

22(®) 

12-15(1) 

21(8) 

21 

M( 10 -’^ Meyr-l) 

1.5-3(9’10) 

20(®) 

11(1) 

to 

1 

00 

to 

1 

bo 

r (pc) 

0.2(11) 

0.2(12) 

0.2(1) 

0.2(18) 

0.3 

ne([S II ]) (cm“®) 

~ioooo(i'‘’*) 

500-1800(1®) 

RilOOO(l) 

^500(18) 

:5i600 

TVfring 

0.06(1®) 

o.i(i®) 

0.5(1) 

0.5-1.0(8) 

1 

2 ;sini( kms“^) 

- 

53 

44 

34 

90 

i{°) 

43 (iR 

65(12) 

60(18) 

50(18) 

60 

Vexpi kms“l) 

10(19) 

30(15) 

2o(20) 

19(18) 

8 

tkin (10"* yr) 

2 

0.7 

1 

1 

3.7 

Mkin/M 

10-20 

7 

45 

250 

40-60 

log(N/H)-tl2 

8.44(1®) 

8.9l(®) 

8.42(1) 

7.51(18) 

8.21 


References: (1) Walborn et al. 1989; (2) Moffat 1983; (3) Walborn & Fitzpatrick 2000; (4) Taylor et al. 2014; (5) Arnett et al. 1989; (6) 
Hendry et al. 2008; (7) Nota et al. 1996; (8) Smith et al. 2013; (9) Blondin & Lundqvist 1993; (10) Martin & Arnett 1995; (11) Panagia 
et al. 1991; (12) Brandner et al. 1997a; (13) Smith et al. 2007; (14) Plait et al. 1995; (15) Brandner et al. 1997b; (16) Mattila et al. 
2010; (17) Jakobsen et al. 1991; (18) O’Hara 2003; (19) Crotts & Heathcote 1991; (20) Hutsemekers et al. 1994. 

Note: ‘Based on the fading of the [O III] A5007 emission line. 


Table 8. Peculiar transverse (in Galactic coordinates) and radial velocities, and the total space velocity of MN18 for two adopted proper 
motion measurements and three distances (see the text for details). 


d 

(kpc) 

fla COS S 
(mas yr“^) 

^5 

(mas yr~^) 

Vl 

(kms“^) 

Vb 

( kms“^) 

Vt 

( kms“^) 

V* 

(kms“^) 

4.4 

-10.3±2.0“ 

-7.6±2.0“ 

-144.3±42.1 

-12.6±41.9 

9.6±1.3 

145.2±42.0 

4.4 

-12.3±2.5® 

-10.3±2.5® 

-209.7±52.1 

-39.5±52.1 

9.6±1.3 

213.6±52.0 

5.6 

-10.3±2.0“ 

-7.6±2.0“ 

-164.4±53.6 

-18.2±53.4 

22.1±1.3 

166.9±53.1 

5.6 

-12.3±2.5® 

-10.3±2.5® 

-247.8±66.4 

-52.4±66.4 

22.1±1.3 

254.2±66.2 

7.1 

-10.3±2.0“ 

-7.6±2.0“ 

-184.3±68.0 

-25.3±67.6 

20.3±1.3 

187.1±67.6 

7.1 

-12.3±2.5® 

-10.3±2.5® 

-290.0±84.1 

-68.7±84.1 

20.3±1.3 

298.7±83.9 


Note: “SPM4.0 (Girard et al. 2011); ^UGAC4 (Zacharias et al. 2013). 


to the stellar space velocity would affect the shape of the 
bipolar nebula, sweeping back the polar lobes, which is def¬ 
initely not the case for the nebula around MN18. One can¬ 
not therefore exclude the possibility that the proper motion 
measurements for MN18 are erroneous (e.g. because of the 
presence of the equatorial ring) and that this star was actu¬ 
ally born not far from its current position on the sky. There 
is some evidence in support of this possibility. Before dis¬ 
cussing this evidence, we note that the small peculiar radial 
velocity of MN18 implies that the distance to this star did 
not change appreciably during its lifetime. Since the pecu¬ 
liar transverse velocity of MN18 is unrealistically high, and 
massive stars did not form in the interarm regions, we con¬ 
clude that the most likely birthplace of MN18 is located in 
the Crux-Scutum Arm, i.e. at d^4 kpc. 

In Section ISTII we noted that the RSS slit crosses two 
other stars on both sides of MN18 (see also Fig. ll3ll . The nor¬ 
malized spectrum of the brightest of these stars, located at 
RA(2000)=15‘'16™42:97, Dec.(2000)=-58°22'19'.'7 and in¬ 
dicated in Fig. [13] by a diamond, is shown in Fig. ll4l In 
the following, we will use for this star its 2MASS name - 
2MASS J15164297-5822197, or J1516-5822, in short. Like 


in the case of MN18, the spectrum of J1516—5822 is dom¬ 
inated by absorption lines of H and Hel. But additionally, 
it also shows weak absorption lines of Hell at AA4686, 5412 
and 6527, which implies that J1516—5822 is of a spectral 
type earlier than B1 (Walborn & Fitzpatrick 1990), i.e. hot¬ 
ter than MN18. The Mgll A4481 line is weak, which is typ¬ 
ical of late O/BO stars. The low quality of the blue part 
of the spectrum precludes us from using traditional clas¬ 
sification criteria. Instead, we used the EW(H 7 )-absolute 
magnitude calibration by Balona & Crampton (1974), which 
for EW(H7)=3.14±0.20 A implies that the spectral type of 
J1516—5822 is either 09 V, BO-1 III or B5-7Ib. The B5“ 
7Ib classification could be rejected because of the presence 
of the He ll lines in the spectrum. Thus, we tentatively clas¬ 
sify J1516—5822 as a O9V/B0-1HI star. EWs, FWHMs 
and RVs of main lines in the spectrum of J1516—5822 are 
summarized in Table[9| 

Detection of a massive and therefore very rare star 
at only «17" from MN18 suggests that both stars could 
be members of the same star cluster. Assuming that 
J1516—5822 is an 09V star and using its 2MASS magni¬ 
tudes (J,H,Ks)=(11.35,10.85,10.63) (Cutri et al. 2003) and 
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Figure 13. Spitzer MIPS 24(left) and SHS Ha+[Nll] images of a field around MN18 with positions of the newly-identified late 
0-/early B-type star 2MASS J15164297—5822197 and star 11 from Lynga (1964a) indicated by a diamond and circle, respectively. The 
orientation of the RSS slit (PA=72°) is shown by a (blue) rectangle. See text for details. 


the photometric calibration of optical and infrared mag¬ 
nitudes for Galactic O stars by Martins & Plez (2006), 
we derived a A'-band extinction towards J1516—5822 of 
Ak = 0.59 mag and a distance modulus of 13.36 mag, 
which corresponds to the distance of 4.7 kpc. The dis¬ 
tance to J1516—5822 would be of «7 kpc if the spectral 
type of this star is BO-1 III. These estimates are consis¬ 
tent with the distance estimate of 5.6ti['2 for MN18 (see 
Section|33J j which supports the supposition that both stars 
might be members of the same cluster. Note that the mean 
RV of lines in the spectrum of J1516—5822 of 145 kms“^ 
is more than twice higher than Wsys of MN18. Although 
this difference might indicate that the two stars are un¬ 
related to each other and are simply projected by chance 
along the same line-of-sight, it could also be understood 
if J1516—5822 is a close massive binary. Note also that 
the SPM4.0 and UCAC4 proper motion measurements for 
this star (pia cos 5 = —28.42 ± 2.89, pis = 45.47 ± 2.92 and 
pia cos 5 = —36.0 ± 3.5, fis = 56.6 ± 3.4, respectively) differ 
significantly from those of MN18. These measurements, how¬ 
ever, seem to be erroneous because they imply a huge space 
velocity of '^1000-2000 kms“^ for any reasonable distance 
to the star. To conclude, the existing data did not allow us 
to prove or disprove the physical association between MN18 
and J1516-5822. 

Finally, we discuss the bright arc-like feature attached 
to the southwestern edge of the southeast lobe of the bipo¬ 
lar nebula around MN18 (see Section[2] and FigsfT] and 1131) . 
Fig.[T3] shows that there is a star within the arc and near 
its apex. This star, located at RA(J2000) = 15*'16"'40\76, 
Dec.(J2000) =—58°23'38”7, is listed in tableV of Lynga 
(1964a) as a member of Lynga 3 under the number 11 
(Lynga3-11 hereafter). We speculate that the arc might be 
created because of interaction between the wind of Lynga 3- 
11 and the nebula (cf. Mauerhan et al. 2010; Burgemeis- 
ter et al. 2013). The enhanced brightness of the arc at the 
place of possible contact with the nebula supports this pos- 


Table 9. EWs, FWHMs and RVs of main lines in the spectrum 
of 2MASS J15164297-5822197. 


Ao(A) Ion 

EW(A) 

(A) 

FWHM(A) 

(A) 

RV 

(km s~^) 

4340 H 7 

3.14±0.20 

15.03±2.53 

-96±16 

4471 He I 

1.00±0.04 

9.58±0.40 

-158±9 

4922 He i 

0.25±0.04 

4.50±0.51 

-163±8 

5876 He i 

0.87±0.04 

6.73±0.27 

-141±6 

6563 Ha 

1.74±0.05 

10.62±0.24 

-146±6 

6678 He i 

0.45±0.03 

7.03±0.17 

-150±6 

7065 He i 

0.71±0.03 

11.16±0.29 

-161±5 


sibility. If the interaction indeed takes place, then Lynga 3- 
11 should be a source of strong wind, i.e. an OB star. 
In the absence of spectroscopic observations of Lynga 3-11, 
we estimate its spectral type using its 2MASS magnitudes 
(J,H,Ks)=(11.07,10.59,10.34) and a photometric calibration 
for Galactic O stars by Martins & Plez (2006), and assum¬ 
ing that this star is located at the same distance as MN18 
(i.e. at d = 5.67)'2 kpc). In doing so, we found that Lynga3- 
11 should be of spectral type of 09.5 V (d=4.4 kpc), 06 V 
(d=5.6 kpc) or 04V/O7.5III (d=7.1 kpc). A spectroscopic 
study of this and other stars listed in Lynga (1964a) is nec¬ 
essary to check whether or not they form the parent cluster 
of MN18. 


7 SUMMARY 

In this Paper, we presented the results of optical spectro¬ 
scopic observations of the blue supergiant MN18 and the 
equatorial ring of its hourglass-like nebula. MN18 was clas¬ 
sified as Blla and its spectrum was modelled using the ra¬ 
diative transfer code cmfgen. We found the stellar effec¬ 
tive temperature of 21.1±1.0 kK and the surface gravity of 
logg=2.75±0.20. Adopting an absolute visual magnitude of 
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Figure 14. Normalized spectrum of 2MASS J15164297—5822197 obtained with the SALT on 2012 March 2 with the principal lines and 
most prominent DIBs indicated. 


Mv^=—6.8±0.5 (typical of B1 supergiants), MN18 has a lu¬ 
minosity of log L/ Lq«5.42±0.3 and resides at a distance of 
«5.6l]^’2 kpc. It was shown that the nitrogen abundance of 
the stellar wind is enhanced by «4 times due to the CNO cy¬ 
cle and that the rotational velocity of MN18 is «100 kms“^. 
The Ha line in the spectrum of MN18 changes from emission 
to absorption on a yearly time scale and shows a noticeable 
variability on a daily time scale. Correspondingly, the mass- 
loss rate estimates based on synthesizing the Ha line show 
a range of values from «2.8 to 4.5 x 10“^ MQyr“^. Our 
photometric observations showed that the B, V and B mag¬ 
nitudes of MN18 changed by less than RiO.l mag during the 
last five years, while its B — H and B — B colours remained 
almost constant. 

Analysis of the nebular spectrum revealed a factor of r; 2 
overabundance of nitrogen, which implies that the equatorial 
ring is composed of processed material. The mass of the ring 
was found to be of Mq, and its systemic and expansion 
velocities were estimated to be of —66.Oil.3 kms“^ and 
7.Oil.2 kms“^, respectively. The latter estimate implies the 
kinematic age of the equatorial ring of Ri3.7 x yr. 

The derived parameters of MN18 were compared with 
those predicted by stellar evolutionary models and a con¬ 
clusion was drawn that this star is in the pre-red supergiant 
phase of its evolution. This conclusion agrees with previ¬ 
ous suggestions that circumstellar nebulae around massive 
stars could be formed already near the end of the main se¬ 
quence phase, and allowed us to suggest that the origin of 
the equatorial ring around MN18 is unlikely to be explained 


within the framework of single star evolution. Based on the 
observed high percentage of binaries among massive stars 
and the possibility that massive stars did not born as rapid 
rotators, but became such because of mass transfer or merg¬ 
ers in binary systems, we suggested that these two binary 
interaction processes might be the main causes for the ori¬ 
gin of nearly critically rotating stars and equatorial outflows 
from them. 

The parameters of MN18 and its equatorial ring were 
also compared with those of four other blue supergiants 
with hourglass-like nebulae in the Milky Way (Sher25, 
HD 168625, [SBW2007] 1) and the Large Magellanic Cloud 
(Sk—69°202). In particular, the comparison revealed that 
the luminosities of these stars span a wide range, from 
log(L/LQ)= 4.7 to 5.8, which implies that the stellar lu¬ 
minosity is not a decisive factor in formation of equatorial 
rings. Also, the nitrogen abundances of the equatorial rings 
around MN18, HD 168625 and [SBW2007] 1 are consistent 
with the possibility that these rings were formed at the end 
of the main sequence stage or soon after it. 

We found that MN18 is among a group of 25 stars sug¬ 
gested in the literature as being members of the star cluster 
Lynga3. The reported age estimates for the cluster, how¬ 
ever, about two orders of magnitude larger than the age of 
MN18, which raises doubts in their association. On the other 
hand, our spectroscopic observations led to a serendipitous 
discovery of a late 0-/early B-type star at an angular sepa¬ 
ration of only «17 arcsec from MN18. The close proximity 
of the two massive and therefore very rare stars suggests 
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that both could be members of the same young star clus¬ 
ter. The distance estimates for these two stars support this 
suggestion. Moreover, we found an indication that another 
member of the putative cluster Lynga 3 is a massive star as 
well, which also points to the possibility that Lynga 3 is a 
young star cluster. Spectroscopic follow-ups of this and other 
stars around MN18 are necessary to prove whether Lynga 3 
is the parent cluster to MN18. 
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